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Abstract: In this note we try to prove the absence of the ghosts in case of the general 
non-linear massive gravity action in Stuckelberg formalism. We argue that in order to 
find the explicit form of the Hamiltonian it is natural to start with the general non-linear 
massive gravity action found in arXiv:1106.3344 [hep-th]. We perform the complete Hamil- 
tonian analysis of the Stuckelberg form of the minimal the non-linear gravity action in this 
formulation and show that the constraint structure is so rich that it is possible to eliminate 
non-physical modes. Then we extend this analysis to the case of the general non-linear 
massive gravity action. We find the corresponding Hamiltonian and collection of the pri- 
mary constraints. Unfortunately we are not able to finish the complete analysis of the 
stability of all constraints due to the complex form of one primary constraint so that we 
are not able to determine the conditions under which given constraint is preserved during 
the time evolution of the system. 
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1. Introduction and Summary 

The first formulation of the massive gravity was performed by Q at least in its linearized 
level as the propagation of the massive graviton above the flat background 1 . Even if the 
theory seems to be well defined at the linearized level there is a Boulware-Deser ghost Q in 
the naive non-linear extension of the Fierz-Pauli formulation. On the other hand recently 
there has been much progress in the non-linear formulation of the massive gravity without 
the Boulware-Deser ghost ||, || and also ]3^, [3c| 2 . 

The Hamiltonian treatment of non-linear massive gravity theory was performed in 
many papers with emphasis on the general proof of the absence of the ghosts in given 
theory. The first attempt for the analysis of the constraint structure of the non-linear 
massive gravity was performed in |^]. However it turned out that this analysis was not 
complete and the wrong conclusion was reached as was then shown in the fundamental 
paper pOf where the complete Hamiltonian analysis of the gauge fixed form of the general 
non-linear massive gravity was performed. The fundamental result of given paper is the 
proof of the existence of two additional constraints in the theory which are crucial for the 
elimination of non physical modes and hence for the consistency of the non-linear massive 
gravity at least at the classical level. 

Then the Hamiltonian analysis of the non-linear massive gravity in the Stiickelberg 
formulation was performed in |^8|, [49t| . Unfortunately the wrong conclusion was again 
reached in the first versions of given papers as was then shown in [^(J where the absence 
of the ghosts in the minimal version of non-linear massive gravity was proven for the first 
time. Then an independent proof of the absence of the ghosts in the minimal version of 
non-linear massive gravity in the Stiickelberg formulation was presented in plf . 

However the proof of the absence of the ghosts in the general form of the non-linear 
massive gravity in the Stiickelberg formulation is still lacking. The difficulty with the 
possible Hamiltonian formulation of given theory is that the action depends on the kinetic 



1 For recent review and extensive list of references, see 

2 For related works, seeJH R |, |, 0, 0, [T| || B ElU 111' !!> H H H H H II H 
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terms of the Stiickelberg fields in a highly non-linear way so that it seems to be impossible 
to find an explicit relation between canonical conjugate momenta and the time derivatives 
of the Stiickelberg fields. On the other hand there exists the formulation of the non-linear 
massive gravity action with the linear dependence on the kinetic term for the Stiickelberg 
fields. This is the form of the non-linear massive gravity action that arises from the original 
one when the redefinition of the shift functions is performed |M 37, 41, 42]. The goal of 



this paper is to perform the Hamiltonian analysis of given action and try to identify all 
constraints. 

As we argued above the main advantage of the formulation of the non-linear massive 
gravity with redefined shift function is that the kinetic term for the Stiickelberg fields 
appears linearly and hence it is easy to find the corresponding Hamiltonian even for the 
general form of the non-linear massive gravity action. Further it is possible to identify 
four primary constraints of the theory where the three ones are parts of the generators of 
the spatial diffeomorphism. Note that the presence of the diffeomorphism constraints is 
the reflection of the fact that we have manifestly diffeomorphism invariant theory. On the 
other hand the fourth primary constraint could be responsible for the elimination of the 
additional non physical mode. However this claim is only true when the requirement of 
the preservation of given constraint during the time development of the system generates 
another additional constraint. Unfortunately we find that the original primary constraint 
cannot provide such additional constraint due to the fact that the Poisson bracket between 
the primary constraints defined at different space points is non zero. For that reason we 
should find another constraint that obeys the property that Poisson bracket between these 
constraints defined at different space points is zero. We find such a constraint in the case 
of the minimal non-linear massive gravity action and we show that this constraint has the 
same form as the primary constraint found in [51|. Then we will be able to show that the 
requirement of the preservation of given constraint during the time evolution of the system 
implies the additional constraint and these two constraints together allow to eliminate two 
non-physical modes. This result agrees with previous two independent analysis performed 
in §(J and in Q. 

Unfortunately we are not able to reach the main goal of this paper which is the proof of 
the absence of the ghosts for the general non-linear massive gravity theory in Stiickelberg 
formalism. The reason is that we are not able to find the primary constraint that has 
vanishing Poisson bracket between these constraints defined at different space points and 
that has the Poisson brackets with another constraints that vanish on the constraint surface. 
This is crucial condition for the existence of the additional constraint. It is rather worrying 
that we are not able to finish the Hamiltonian analysis for the general non-linear massive 
gravity action especially in the light of the very nice proof of the absence of the ghosts in 
case of the gauge fixed non-linear massive gravity action [40|. However there is a possibility 
that the proof of the absence of the ghosts for general non-linear massive gravity action in 
Stiickelberg formalism could be found in the very elegant formulation of the massive and 
multi metric theories of gravity presented in [23]. We hope to return to this problem in 
future. 

The structure of this paper is as follows. In the next section (pi) we introduce the 
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non-linear massive gravity action and perform the field redefinition of the shift function. 
We also we find the Hamiltonian formulation of the minimal form of the non-lineal massive 
gravity in Stiickelberg formulation with redefined shift functions and we find that given 
theory is free from the ghosts. Then we extend this approach to the case of the general 
non-linear massive gravity theory in the section (??). We identify all primary constraints 
and discuss the difficulties that prevent us to finish the complete Hamiltonian analysis. 



2. Non-linear Massive Gravity with Redefined Shift Functions 

As we stressed in the introduction section the goal of this paper is to perform the Hamil- 
tonian analysis of the general non-linear massive gravity with presence of the Stiickelberg 
fields. It turns out that it is useful to consider this action with redefined shift functions 
1 38, |3^, f|l|, |42[| . More explicitly, let us begin with following general form of the non-linear 
massive gravity action 

r 3 

S = M 2 p d 4 xy^[( 4 ) J R + 2m 2 ^/3 n e n (7r T 7)] , (2-1) 

n=0 

where efc(A) are elementary symmetric polynomials of the eigenvalues of A. For generic 
4x4 matrix they are given by 

e (A) = 1 , 
ei(A) = [A] , 

e 2 (A) = i([A] 2 - [A 2 ]) , 

e 3 (A) = i([A] 3 -3[A][A 2 ]+2[A 3 ]) , 

e 4 (A) = ^([A] 4 -6[A] 2 [A 2 ] + 3[A 2 ] 2 + 8[A][A 3 ]-6[A 4 ]) , 
e fc (A) = , for k > 4 , 

(2.2) 

where A^ is 4 x 4 matrix and where 

[A] = TrA^ . (2.3) 

Of the four (3 n two combinations are related to the mass and the cosmological constant 
while the remaining two combinations are free parameters. If we consider the case when the 
cosmological constant is zero and the parameter m is mass, the four /3 n are parameterized 
in terms of the 03 and 04 as || 

/?n = (-l) n Q(4-n)(3-n)-(4-n)a 3 + a 4 ) . (2.4) 

The minimal action corresponds to fi% = (3% = that implies «3 = 04 = 1 and consequently 
A) = 3 ,pi = -1. 
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We consider the massive gravity with that is manifestly diffeomorphism invariant. This 
can be ensured with the help of 4 scalar fields (j> A , A = 0, 1, 2, 3 so that 

sTUp = gTdv^dptA . (2.5) 

Then we have 

at2 . = [ -/oo + N l fi Q -f j + N l fij \ 

9 1 \ N 2 g il f l0 - iV*(-/oo + N l f lQ ) iVV/ij - N*{-f oj + J ' 1 ' ] 



where we also used 3 + 1 decomposition of the four dimensional metric g^ u [ [15; , 46] 



5oo = ~N + Nig l3 Nj , g i = N { , g^ = , 

~oo = _±_ ~oi = ^_ ~g _ ij jgjg 
9 N 2 ' 9 N 2 ' 9 9 N 2 ' 

(2.7) 

Let us now perform the redefinition of the shift function N l that was introduced in fj7| , |38| , 



H HI 

TV* = Mff + f ik f ok + ND^nP , (2.8) 

where 

x = 1 - n%n J ' , M 2 = -/oo + / fc/ H //o (2.9) 
and where we defined f tJ as the inverse to /y in the sense 3 

= 5/ . (2.10) 

Finally note that the matrix D l ■ obeys the equation 



kj 



(2.11) 



and also following important identity 



fikD j = fjkb\ . 

(2.12) 

Let us now concentrate on the minimal form of the non- linear massive gravity action. 
Using the redefinition (p.8| ) we find that it takes the form 

S = Ml J cPmttlNJgKij&MKu + Ny/gR - JgMU - 2m 2 (N^^D\ - 3N^g)] , 

(2.13) 

where 

U = 2m 2 VE , (2.14) 



s Note that in our convention f ik coincides with ( 3 / -1 )^ presented in Q § |1[ |§. 
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and where we used the 3 + 1 decomposition of the four dimensional scalar curvature 



^R = K ij g» kl K kl + R, (2.15) 
where R is three dimensional scalar curvature and where 

gijkl = itfkgjl + gilgjkj _ g ij g kl (2 



with inverse 



' ' — -iklmn 1 



Gijki = -^{gikgji + gugjk) - ^gijgki , QijkiQ mn = -^t^ + 8? 8™) • (2.17) 

Note that in ( p,15|) we ignored the terms containing total derivatives. Finally note that 
Kij is defined as 

Kij = ^(dtgij - ViNj{n,g) - VjN t (h,g)) , (2.18) 

where iVj depends on h 1 and g through the relation ( |2.8| ). 

At this point we should stress the reason why we consider the non-linear massive gravity 
action in the form (l2~Ig| ). The reason is that we want to perform the Hamiltonian analysis 
for the general non-linear massive gravity action written in the Stiickelberg formalism. It 
turns out that the action ( p. 13 ) has formally the same form as in case of the general non- 
linear massive gravity action when we replace U with more general form whose explicit 
form was determined in [[37], [38], ^2j. On the other hand the main advantage of the 
action ( 2.13j ) is that it depends on the time derivatives of (j) through the term M and 
that this term appears linearly in the action ( 2.13 ). We should compare this fact with the 
original form of the non-linear massive gravity action where the dependence on the time 
derivatives of <fi A is highly non-linear and hence it is very difficult to find corresponding 
Hamiltonian. 

Explicitly, from (|1|) we find the momenta conjugate to N, h % and gij 

n N « ,7Ti w ,7r« = M 2 p ^ kl K kl (2.19) 
and the momentum conjugate to 4> A 

(2.20) 

where 

Ki = -2g lk V^ . (2.21) 
It turns out that it is useful to write M 2 in the form 



M 2 = -d t( t> A M AB d t ct> B , M AB = VAB - difaFdjfo 



(2.22) 
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where by definition the matrix Mab obeys following relations 

MabV BC M C d = Mad , detM A B = 1 



(2.23) 



together with 

d l <P A M A B = di<t> B - d^d^Af'd^B = . (2.24) 
With the help of these results we find 

Pa + KipdrfA = (h l TZi + M^U)-^M A Bd t (j) B (2.25) 

and consequently 

M 2 

M 2 = -d t <t> A M A Bd t <t> B = - {hiUi + M2 )2 (pa + Kifidj^n^ipB + Kifid^B) 

(2.26) 

which however implies following primary constraint 

S p = [fjKi + M 2 p ^Uf + (p A + Kif ij d j( t>A)(p A + Kif j d j( t> A ) « . (2.27) 
Note that using ( |2.24| ) we obtain another set of the primary constraints 

d^ A u A = di4> A VA + Hi = « . 

(2.28) 

Observe that using ( |2.28 ) we can write 

p A + TZifd^A = MacV CB Pb + >:j iJ <)jO A (2.29) 

so that we can rewrite E p into the form 

X p = (h i K i + M*^gU) 2 + Pa M Ab pb + H% , (2.30) 

where H l are functions of the phase space variables. As a result we see that it is natural 
to consider following independent constraint E p 

S p = (n^ + M 2 v ^Xjf + p A M AB p B ~ . (2.31) 

We return to the analysis of the constraint S p below. 

Now we are ready to write the extended Hamiltonian which includes all the primary 
constraints 



H E = j d 3 x(NC + v n tt n + v% + fi p S p + Q%) , (2.32) 



where 



Cq = -^TT^g i:ikl 7T kl - Mp^/gR + 2m 2 M 2 ^ViD i i - 6m 2 M 2 ^ + D^nPTli 
v Q p 



(2.33) 
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and where we introduced the constraints Si defined as 

Si = Si + difj^i + dj {h j 7Ti) . (2.34) 

Note that Si is defined as linear combination of the constraints Si ~ together with the 
constraints 7Tj ~ 0. 

To proceed further we have to check the stability of all constraints. To do this we have 
to calculate the Poisson brackets between all constraints and the Hamiltonian He- Note 
that we have following set of the canonical variables gij,TT 1 ^ ,(f) A ,pA-,n l ,TTi and N,ttn with 
non-zero Poisson brackets 

{^•(x),7r M (y)} = I(^ + ^)5(x-y) , {^(x),p B (y)} = $tf(x - y) , 
{iV(x),^v(y)} = ff(x-y) , {^(x), ^(y)} = <5^(x - y) . 

(2.35) 

Now we show that the smeared form of the constraint Si ~ 

T S (C) = [ d 3 x(% (2.36) 



is the generator of the spatial diffeomorphism. First of all using ( |2.35| ) we find 

{T s (C),n k } = -Cd i n k +^d j C k 

which is the correct transformation rule for n\ Then using fl2.35| ) we find 

{TsiN*),^} = -d l N i U j - S'OjUj ■Rjitj.X 1 , 
{T s (N l ), PA } = -N%p A - m'pA , 
{T s {N l ),(p A } = -N l d^ A , 

{T s (N% 9ij } = -N k d k g ij -d i N k g kj -g lk d j N k , 
{TstN*),^} = -d k (N k ir ij ) + d k N l Tr kj + ir tk d k N j , 
{T s (N l )J tj } = -N k d k fij - diN k f kj - f lk djN k , 
{T s (N 1 ),tt' 1 } = -3*iVV - .V • fyiVV , 



(2.37) 



(2.38) 



that are the correct transformation rules of the canonical variables under spatial diffeo- 
morphism. To proceed further we need the Poisson bracket between Ts(iV 4 ) and D l -. It 
turns out that it is convenient to know the explicit form of the matrix D % ■ |37], |38], 41, 42] 



D\ = ^Jg^fmnQMQ- 1 )^ , (2.39) 

where 



Q l j = xPj + n*n%- , (Q- 1 )^ = ±(6* - hPh m f mq ) . (2.40) 
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Using this expression we can easily determine the Poisson brackets between Ts(iV l ) and 
D k j. In fact, by definition we have 

{t^),*?*,} = -N m d m Q\ + d m N k Q"t - Q k n d t N n , 
{T 5 (iV i ),g l J Q i fc } = -N m d m {Q^Q\) + d m N*Q l J Q\-Q i 3 Qi m d k N m . 

(2.41) 

Using ( |2.39| ) and the results derived above we find 

{TsCJV*),^-} = -N m d m D i J + dmNW™ - D^djN™ . (2.42) 

Collecting all these results and after some calculations we find 

{TsCJV^Cb} = -N m d m c - d m N m C , 
{T 5 (iV l '),S p } = -N m d m Z p - d m N m Y> p . 

(2.43) 

Note also that it is easy to show that following Poisson bracket holds 

{TsiN*), T s (M j )} = T s (N j d j M i - M^djN 1 ) . (2.44) 

Now we are ready to analyze the stability of all primary constraints. As usual the require- 
ment of the preservation of the constraint tt^ ~ implies an existence of the secondary 
constraint Cq ~ 0. However the fact that Cq is the constraint immediately implies that 
the constraint E$ ~ is preserved during the time evolution of the system, using ( 2.45| ) 
and (2.44). Now we analyze the requirement of the preservation of the constraints 7Tj ~ 
during the time evolution of the system 

dm = {m, H E } = - (n p 5f + d{D g^ J j {n k - 2m 2 M 2 p ^f km n m ^ = . (2.45) 



It turns out that the following matrix 



n ptf + Ji = ( 2 - 46 ) 



cannot be solved for Q p and hence we have to demand the existence of following secondary 
constraints H, EE Hi 

Ci = TZi jP^fijfiP « . (2.47) 

yx 

Finally we have to proceed to the analysis of the time development of the constraint 
T, p « 0. However it turns out that it is very difficult to perform this analysis for T, p due 
to the presence of the terms that contain the spatial derivatives of <j) A . Then the explicit 
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calculation gives {S p (x), S p (y)} 7^ 0. For that reason we proceed in a different way when 
we try to simplify the constraint S p . Using Cj and Sj we find that 

A = J + F l Si + G J C 4 , 

(2.48) 

where 

,4 = PAdi^f'dj^ps , (2.49) 

and where F l , G % are the phase space functions whose explicit form is not important for us. 
Then with the help of ( p. 48 ) we express n % fijhi as a function of the phase space variables 



p A ,4> A and gij,Tr l 3 

A-F^i- G l C t 
(A - F^i - G'Ci) + 4m 4 M 4 ff ' 

In the same way we obtain 



= ^ T^-^IaVTaZ • ( 2 - 5 °) 



n l ll % = A ° tCl + , (2.51) 

'(A - F^i - G%) + 4m 4 M 4 9 



and 

n l = J ^ PAJ + F l £; + G*Cj , (2.52) 

/ A + 4m 4 M 4 ff 



where again F 1 , G % are phase space functions whose explicit form is not needed for us. 
Now using these results we find that the constraint £ p takes the form 



v _ (A - - G l Q + 4m 4 M 4 g)4m 4 M 4 g A 
^ p " ^-F^-G*C 4 + 4m 4 M 4 5 + ^ ^ +PAP 



+ 4m 4 M 4 9 + f^Si = 4m 4 M 4 5 S p + H% , 



(2.53) 



where we introduced new independent constraint S 



^4 

S„ = - +1 = (2.54) 

p 4m 4 M 4 5 v ; 



that has precisely the same form as in [51]. Note that the constraint £ p has the desired 
property that |s p (x), S p (y)| =0. As a result we see that it is more natural to consider 

S p instead of E p as an independent constraint. Then the total Hamiltonian, where we 
include all constraints, takes the form 

H T = j d 3 x(NC + v N ir N + uV* + Sl p t p + U% + r%) . (2.55) 
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Now we are ready to analyze the stability of all constraints that appear in ( |2.55 ). First of 
all we find that ~ is automatically preserved while the preservation of the constraint 
7Tj gives 



-2m l V — (f l3 - f ik n k j jX n l ) = -A^^P 
vx 



(2.56) 



By definition 



det(/y - fikn k fuh l ) = xdet/y / 



(2.57) 



and hence the matrix A^g. is non-singular. Then the only solution of the equation ( 2.56| ) 
is P = 0. 

As the next step we proceed to the analysis of the stability of the constraint S p . As is 



clear from (2.54) we have 



{s p (x),S p (y)}=0. 
Then the time evolution of given constraint takes the form 



(2.58) 



d t t p = [£ P ,H T } « J d 3 xAT(x){S p ,C (x)} 



(2.59) 



using the fact that E p does not depend on h l together with T l = and also the fact that 
E p is manifestly diffeomorphism invariant. 

In order to explicitly determine ( |2.59 ) we need following expression 



Sfij 



^D\r - -±=h i f lm ~ 



(2.60) 



Then after some calculations we obtain 



S p , / d 6 xNC \ = 2d i [ND l j ]n :) T, p + 



M*m*g 



+ N —D' 



Pa9, 



N 



C k d^ A 



d l [n j PA ]PA , 2m 2 M ? 2 



M£m A g 



p A dj(f) A di 



~Sf< 



N 



ij 



3 2m p m 4l g M p m^g 



+ T7Iz£.D l kPAdi[^Vif kj d^ A ] « iVS 



n 



(2.61) 
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In order to simplify E 7/ further we use ( ET501) , (|2.51|) together with (gjp to make E^ 7 
independent on n\ In fact, using ( 2,39| ) ( pT4Q|) we obtain 

Qp = A + 4 M4 m 4 g (^M*g6™ + B^pjJ^&^PB) , 
A + AM^m 4 g 1 

(2.62) 



up to terms proportional to the constraints Cj, Ej. With the help of these results we obtain 

^ J 2m A m 4 g M A m A g p 

= t 1 / + i^S, + G^Cj , 

(2.63) 

where x and Z)* are functions of pA,dj(f) A and g through the relations ( 2.50^ ,( 2.51 ) and 
(2.52). Now we see from (2.59) that the time evolution of the constraint E p ~ is obeyed 
on condition when either N = or when E 7 / = 0. Note that we should interpreted N as 
the Lagrange multiplier so that it is possible to demand that N = on condition when 
Ti 11 ^ on the whole phase space. Of course such a condition is too strong so that it is 
more natural to demand that S 77 ~ and N ^ 0. In other words E 7 / ~ is the new 
secondary constraint. 

In summary we have following collection of constraints: tin ~ , 7Tj ~ ,Cq ~ 0,Cj ~ 
0, Ej « 0, E p « 0, E 7/ w 0. The dynamics of these constraints is governed by the total 
Hamiltonian 

H T = J d 3 x(NC + v N ir N + v^i + %t p + fi 77 E 7/ + O^E* + Y%) . (2.64) 

As the final step we have to analyze the preservation of all constraints. The case of 7Ttv ~ 
is trivial. For 7Tj ~ we obtain 

d m (x) = K(x),tf T } = j d 3 y(r J (y){^(x),C J (y)} + ^ / (y){vr J (x),S 7/ (y)}) = 
= PA^ Cj .(x) = 

(2.65) 

due to the crucial fact that E 7 / does not depend on n\ This is the main reason why we 
introduced E 7/ instead of E 77 . Then as we argued above the only solution of the equation 
is T l = 0. Now the time development of Cj is given by the equation 

d 3 x (N(y) {C,(x),C (y)} + w*(y) {C,(x), ^(y)} + 
n p (y) jc^E^y)} + fi{ J (y) {^(x), S^(y)}) 

(2.66) 



+ 
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and the time development of the constraint T, p is governed by the equation 
d t t p (x) = {t p (x),H T }^ y ( i 3 x^ / (y){s p (x),^ / (y)} . 

(2.67) 

As follows from the explicit form of the constraint Tip 1 we see that |llp / (x), S p (y)| is 
non-zero and proportional also to the higher order derivatives of the delta functions. As a 
consequence we find that the only solution of the equation above is Vlp 1 = 0. Further we 
analyze the time evolution of the constraint Yip 1 

3 x (iV(y) {^(x),C (y)} + O p (y) {sJ'(x),E p (y)}) = . 

(2.68) 



d 3 ~ 



Now from the last equation we obtain Q p as a function of the phase space variables and 
N, at least in principle. Then inserting this result into the equation for the preservation of 



Ci ( 2.66 ) we determine v 3 as functions of the phase space variables. Finally note also that 
the constraint Co is automatically preserved due to the fact that T l = Q^ 1 = and also the 
fact that {Co(x), Co(y)} ~ as was shown in |40]| . 

In summary we obtain following picture. We have five the first class constraints ~ 
, Co ~ , Sj together with eight the second class constraints 7Tj , Cj and 
E p » , Tp 1 « 0. The constraints 7Tj m together with d ~ can be solved for 7Tj 
and n\ Then the constraint S p can be solved for one of the four momenta pa while the 
constraint Y 11 can be solved for one of the four (f> A . As a result we have 12 gravitational 
degrees of freedom <?jj,7T* J ,6 scalars degrees of freedom together with 4 first class constraints 
Co ~ , Sj ~ 0. Then we find that the number of physical degrees of freedom is 10 which 
is the correct number of physical degrees of freedom of the massive gravity. 

3. General Non-Linear Massive Gravity Action 

Let us try to apply the procedure performed in previous section to the case of the general 
non-linear massive gravity whose action takes the form 

S = Ml j d 3 xdt[N^jKijg ijkl K kl + y/gNR + 2m 2 ^MU + 2m 2 N^V] , (3.1) 

v = /3 + hVsb'i + ^Vf[b i i b i j + b^bW + 



where 



(3.2) 
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Following the analysis performed in the previous section we find the extended Hamiltonian 
in the form 

H E = j d 3 x(NC + v n tt n + uSr< + ft p £ p + Sl%) , (3.3) 

where 

Co = -^^Gijki^ 1 ~ M*y/gR - 2m 2 M 2 p V + KiD 1 ^ , (3.4) 

and where the primary constraint £ p takes the form 

£ p : (ffKi + 2m 2 M 2 ^U) 2 + p A M AB p B ~ . (3.5) 

As the next step we should analyze the stability of all primary constraints. As in previous 
section we find that the stability of the constraint tin ~ implies the secondary constraint 
Co ~ while the stability of the constraints 7Tj implies set of the secondary constraints Cj 

§7|, m 0, n 



d = Ki- 2m 2 ^g^ bi4 + f3 2 Vz(6iD m m - £){)+ 



-fc^/f (h>\{D m m B n n - D m n D n m ) + ~ D\D m ^ 



(3.6) 



Now we come to the key point of the analysis which is the requirement of the preservation 
of the constraint S p ~ during the time evolution of the system. This is very complicated 
expression which depends on the all phase space variables. Remember that in the minimal 
case we expressed n l as functions of gij,PA and 4> A . As a result we found that E p can be 
expressed as a linear combination of Cj, Sj and S p where S p obeys an important property 
|e p (x), S p (y)| = O.It would be certainly nice to repeat the same procedure in the case 
of the constraint S p given in (3.5). Unfortunately we are not able to solve the constraint 



Cj in order to express n l as a function of IZi. Consequently we are not able to express 
the constraint S p as a linear combination of the constraints Cj and possibly Si, Co and the 
new constraint S p where |s p (x), S p (y)| = 0. In other words, despite of the fact we are 
able to find the primary constraint S p ~ we are not able to determine the additional 
secondary constraint which is necessary for the elimination of two non-physical phase space 
modes. As a result the main goal of this paper which was the proof of the absence of the 
ghosts in the general non-linear massive gravity action in Stiickelberg formalism cannot be 
completed. 
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